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Abstract 

We solve the 0{N) vector model at large A on a squashed three-sphere with a 
conformal mass term. Using the Klebanov-Polyakov version of the AdS 4 /CFT 3 corre¬ 
spondence we match various aspects of the strongly coupled theory with the physics 
of the bulk AdS Taub-NUT and AdS Taub-Bolt geometries. Remarkably, we find that 
the field theory reproduces the behaviour of the bulk free energy as a function of the 
squashing parameter. The 0{N) model is realised in a symmetric phase for all finite 
values of the coupling and squashing parameter, including when the boundary scalar 
curvature is negative. 


1 Introduction and summary 

The AdS/CFT correspondence posits the existence of exact dualities between large N field 
theories and string theory on asymptotically Anti-de Sitter (AdS) spacetimes [1, 2, 3]. As is 
well known, these gauge-gravity dualities typically relate two complex theories, each being 
tractable only in some limiting region of parameter space. A somewhat simpler version of 
such a duality was recently proposed by Klebanov and Polyakov [4] wherein the large N field 
theory is exactly solvable. Specifically, the Klebanov-Polyakov correspondence conjectures 
a duality between a theory of massless higher spin gauge fields in AdS 4 spacetime on the 
one hand, and the singlet sector of the critical 0{N) vector model at large N in three 
dimensions on the other. 

Motivated by the Klebanov-Polyakov proposal, in this paper we will solve the 0{N) 
vector model at large N on homogeneous, compact manifolds known as squashed three- 
spheres. The squashed three-spheres are one parameter deformations of the round three- 
sphere which introduce an anisotropy whilst preserving homogeneity. Importantly, they are 
the conformal boundary of certain asymptotically locally Anti-de Sitter spacetimes known 
as AdS Taub-NUT and, for certain values of the squashing parameter, of another geometry 
called AdS Taub-Bolt [5, 6]. Both these backgrounds are Einstein manifolds with constant 
negative curvature. The holographic correspondence then suggests an equivalence between 
higher spin gauge theory on the AdS Taub-NUT/Bolt geometries, and the 0{N) vector 
model at large N on the boundary geometry. 

These backgrounds were studied in the context of the AdS 4 /CFT 3 correspondence in 
[5, 6, 7, 8, 9, 10, 11, 12] . However, in those works little could be said about the dual field 
theory because in the AdS^ x 5^ version of the gauge-gravity correspondence [1] the field 
theory is a strongly coupled conformal field theory which remains largely unknown. This 
is no longer the case for the Klebanov-Polyakov version of the correspondence, where the 
field theory is tractable at all values of the coupling, weak and strong. 

We will study the 0{N) theory on a squashed three-sphere as a function of the squashing 
parameter and a dimensionless coupling constant. In the absence of a complete formulation 
of the theory of higher spin gauge fields on the bulk geometries of interest, we compare exact 
results from field theory with semiclassical properties of Einstein gravity in the bulk [5, 6, 7, 
8, 11]. One of the main results of our study is a remarkable qualitative agreement between 
the free energy of the field theory at strong coupling and that of the bulk gravitational 
theory, as a function of the squashing parameter. This result is encapsulated in figures 3 
and 4 below. 
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The squashed sphere is an bundle over S'^ where the squashing parameter a is related 
to the periodicity of the fibre. The periodicity of the fibre is often interpreted as the 
inverse of a temperature, giving T ~ \foi for large a. This then allows a physical under¬ 
standing of the fact that the bulk gravity theory with a squashed three-sphere as boundary 
exhibits a first order phase transition of Hawking-Page type [13] as this temperature is 
increased beyond a critical value [5, 6]. In particular, for fixed a < acrit the semiclassical 
quantum gravity partition function is dominated by the AdS Taub-NUT geometry. Be¬ 
yond ttcrit however, the partition function is dominated by the AdS Taub-Bolt background. 
The AdS Taub-Bolt geometry should be thought of as a Euclidean black hole with a NUT 
charge. 

Our field theory result for the free energy of the 0{N) vector model at strong coupling 
will agree with the behaviour of AdS Taub-NUT for a < cvcrit and with that of AdS Taub- 
Bolt at large values of a, but with a smooth crossover between the two regimes instead 
of a first order phase transition. We view the qualitative agreement away from a = acrit 
as a positive test of the proposed holographic duality. We will discuss the possibility that 
accounting for the infinite massless higher spin degrees of freedom in the bulk will smooth 
out the first order Hawking-Page transition. 

It is not surprising that the nature of the Hawking-Page transition here is qualitatively 
different to the cases in which it is dual to confinement/deconfinement transitions in large N 
Yang-Mills theory [2, 14, 15, 16]. These transitions are accompanied by a drastic change in 
the number of degrees of freedom of the boundary field theory. The natural phase transition 
one might expect in the 0{N) model is a symmetry breaking transition as one varies a. 
However, we will show explicitly that the large N limit of the 0{N) model on a compact 
space is always realized in an 0(N)-symmetric phase at finite coupling, with any possibility 
of spontaneous symmetry breaking of 0{N) 0{N — 1) precluded by finite volume effects. 

Let us now briefly summarise a few further features and implications thereof of our field 
theoretic results. In flat space, the critical 0{N) vector model in three dimensions exhibits 
a renormalisation group flow from a free fixed point in the ultraviolet (UV) to an interacting 
fixed point in the infrared (IR), see for example [4, 17, 18]. This flow is induced by the 
quartic interaction which is a relevant operator with a coupling constant that has mass 
dimension one. A compact space such as the squashed sphere introduces a second scale, 
thus providing a dimensionless coupling constant. The large N limit leads to an exactly 
solvable theory for any value of this dimensionless coupling constant. 

We consider the theory with a conformal mass term so that in the ultraviolet the theory 
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approaches a free conformal fixed point, presumably corresponding to the asymptotically 
AdS region of the bulk. Interestingly, the conformal mass term becomes negative for a range 
of values of the squashing, as the scalar curvature of the squashed sphere becomes negative. 
One of our results will be that nonetheless the theory is always realised in a symmetric 
phase for all finite values of the coupling and squashing. 

In addition to matching the behaviour of the free energy with squashing as we have 
described above, we make two further connections between the physics of the 0{N) model 
and the bulk theory 

• A negative curvature of the boundary is often associated with instabilities in the 
bulk geometry, as reviewed in [11]. Our field theory results suggest that there is no 
instability of the classical bulk AdS Taub-NUT geometry in the higher spin gauge 
theory when the scalar curvature of the conformal boundary becomes negative. We 
will perform a check of this statement by considering fluctuations of a conformally 
coupled scalar field about the bulk background. 

• We calculate the condensate (<1>-<I>) in the 0{N) model. This vacuum expectation 
value is dual, in the Klebanov-Polyakov correspondence, to the normalisable mode of 
a conformally coupled scalar field in the bulk, ip. This scalar field is not turned on 
in the AdS Taub-NUT or AdS Taub-Bolt geometries. Consistent with this fact we 
find (<h-$) ^ 0 at strong coupling. For boundaries other than the round sphere, this 
condensate does not vanish at finite coupling, predicting curvature corrections to the 
bulk AdS Taub-NUT/Bolt geometries, whilst AdS itself is protected from corrections. 

Given that the nonperturbative large N resummation can be implemented in field theory 
for arbitrary values of the dimensionless coupling, we can in fact analyse the theory both 
at arbitrarily weak coupling and at infinitely strong coupling. In these two limits we find 
analytically the following field theory physics 

• The theory at arbitrarily weak dimensionless coupling, aX 0^, exhibits a phase 
transition at a value of the squashing parameter where the boundary curvature turns 
negative. The order parameter for this transition is the condensate A(<I>-<h) which 
vanishes for negative curvature while remaining nonzero and of order N at positive 
curvatures. Interestingly, for any finite value of the coupling this phase transition gets 
completely smoothed out, as we stated above. Such nonperturbative effects may be 
relevant more generally for attempts to extrapolate perturbative phase transitions to 
strong coupling [14, 15]. 
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• At large squashing parameter the free energy of the AdS Taub-Bolt geometries is 
known to depend linearly on the squashing parameter [5, 6]. We find that the 0{N) 
model exhibits the same linear dependence at both strong and weak coupling. The 
ratio of the free energy at large squashing at weak and strong coupling is 4/5, rather 
analogously to the well-known 3/4 factor [19] for AA = 4 super Yang Mills, except that 
here we can calculate both limits in field theory. 

The physics of interacting quantum fields on a squashed three-sphere received attention 
at various points in the eighties following the growth of interest in quantum fields on curved 
backgrounds. These works included considerations of quantum induced symmetry break¬ 
ing through calculation of the one loop effective potential [20, 21, 22, 23, 24]. However, 
although one loop effective potentials can be useful, they are insufficient to understand 
the phase structure of field theories. Firstly because when quantum corrections become of 
the same order as the classical mass term then typically perturbation theory is not reli¬ 
able. Secondly, the one loop minima of the effective potential are generically located at a 
nonperturbative mass scale, again invalidating the reliability of the perturbative analysis. 
Our results described above for negative values of the conformal mass term highlight the 
inadequacy of a perturbative treatment. 

The layout of this paper is as follows. Section 2 introduces the 0{N) model on a general 
background, including effects due to compactness of the spacetime. Section 3 specialises to 
the squashed three-sphere and studies the conformally coupled theory. Section 4 connects 
the field theory results with the dual gravitational description using the Klebanov-Polyakov 
correspondence. The appendices contain technical calculations of the zeta function on a 
squashed three-sphere, as well as a check on our calculations. 


2 The 0{N) model at large N 


We begin by reviewing known results on the interacting 0{N) scalar field theory on a D 
dimensional spacetime A4 (see e.g. [17]). We also note some subtleties arising from the 
fact that the spacetime A4 is compact. The large N limit of the field theory will allow 
us to analyse the exact effective action for this theory and thus draw reliable conclusions 
about the dynamics of the quantum theory. In Euclidean signature the 0{N) model in D 
dimensions has the classical action 





( 1 ) 
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where is an iV-component field transforming as a vector under 0{N) rotations. The 
interacting scalar field theory is super-renormalisable in dimensions less than four. The 
explicit N dependence in (1) is necessary for the theory to posses a well defined large N 
limit. In the action (1) the mass term includes a possible coupling to the scalar curvature 
of the background. 

To solve the large N theory one introduces an auxiliary scalar field a which permits a 
rewriting of the action as 


/■ 


= I 


XN 


-Vd>-V^> + - (m^ + Ait) $ - —a^ 

9 9 ' ' A 


( 2 ) 


The field a is the composite operator, 


Na = $ • $ . 


(3) 


To derive the large N effective potential we introduce a homogeneous background expec¬ 
tation value for the 0{N) held and integrate out all inhomogeneous fluctuations about this 
conhguration. Without loss of generality we may always choose to rotate this expectation 
value into the top component of the 0{N) vector, so that 

^ = (\/iV0 + (50, VTI, 7r2,..., ttat-i) , (4) 


where 0 is a homogeneous background and (50 and if are the quantum huctuations. As 
we will see shortly, the explicit factor of ^/N in (4) is the correct scaling behaviour for 
VEVs in the large N theory so that 0 ~ 0{N^). At this point it is worth noting that the 
backgrounds we consider in this paper, squashed three-spheres, are homogeneous and hence 
the notion of an effective potential for homogeneous fields makes sense. 

Introducing a homogeneous background might appear to break the 0{N) symmetry to 
0{N—1). However, the path integral must include an integral over the vacuum manifold and 
we will see that this implies that symmetry breaking does not actually occur on a compact 
space. In fixing the 0{N) symmetry for the homogeneous modes in (4) we have effectively 
transformed into polar coordinates for $ with \/iV0 playing the role of the radial coordinate. 
Upon integrating over the angular directions in field space the partition function picks up 
an extra factor given by the volume of the vacuum manifold 0{N)/0{N — 1) = 

T[{N - I)/2] = exp {N/2 [l + Invr + ln(0VM)]) , (5) 

where fx is an arbitrary scale introduced to keep the arguments of the logarithms dimen¬ 
sionless. We note that there is a curious cancellation between the term coming from 
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the radius of the A^-sphere and from the large N limit of the Gamma function. The result 
of this cancellation is that the extra contribution to the partition function (5) has the same 
N dependence as the remaining terms in the effective action in the large N limit. 

In the large N limit the effective potential for the homogeneous fields (/> and a, after 
integrating out tt and including the contribution (5), is ^ 


= -{m +Xa}(p - -a + 


1 


N 


2\olM 


In det ' 


—□ + + Act 

^2 


1 


2VolM 


1 + In TT + In ■ 


J 


( 6 ) 


where VolA4 is the volume of the spacetime and ^ is again an arbitrary dimensionful scale 
which should be interpreted as a sliding renormalisation scale. The prime in det^ denotes 
the fact that the integration over the tt fluctuations did not include the constant modes, 
which we have dealt with separately. 

In the large N limit only the saddle point configuration obtained by extremising (6) 
contributes to the partition function. Extremising Ves{(j), a) with respect to the background 
fields (p and a we then hnd the vacuum equations 


+ Act) 


1 

VolM ’ 


(7) 


and 


{rP^ - ^) + 


1 


VolM 


:TT' 


1 


—□ + m2 + Act 


= 0 . 


( 8 ) 


The second of these equations is in fact, in the large N limit, just the vacuum expectation 
value of the operator equation (3) which would yield a = + (7f2)/iV. Here vf refers to 

the {N — 1) “pions” of equation (4). Thus the functional trace computes the renormalised 
vacuum expectation value of the composite operator vf • tt, encoding the contribution of the 
quantum fluctuations about the ground state. 

At this stage it is convenient to define the “effective pion mass” 


m^ = m^ + Xa , 


(9) 


which is precisely the mass of the {N — 1) “pion” fluctuations in (4). From equation (8) 
we see that this effective mass incorporates a self-consistent resummation of the quantum 
fluctuations about the vacuum state. This is a consequence of the large N limit which 
resums the so-called cactus diagrams of the theory. 

^The contribution from integrating out the fluctuations Sip is subleading in the 1/A^ expansion. 


6 



The effective mass depends on whether the theory is realised in a symmetric phase 
or not. In an O(A^)-symmetric phase is nonvanishing. In a symmetry broken phase we 
would expect to vanish, thus implementing Goldstone’s theorem. That is, in a symmetry 
broken phase there would be — 1 massless Goldstone bosons. In equation (7) we see that 
cannot be zero. To achieve this we would need to take VolAd or to infinity. We will 
see below that taking cfp' to infinity requires taking the coupling to zero. Thus we find that 
symmetry breaking cannot occur for the 0{N) model in a compact spacetime at large N 
and finite coupling. 

Equations (7) and (8) may be rewritten as a gap equation for 

( 10 ) 

Note that the trace now includes the constant mode, as will all remaining traces in this 
paper. Once we have solved (10) for we may evaluate the effective potential at the 
extremum to find the action 

= ^ (m^ - mlf -Hn det f ^ ln(;U^Vol>l) -hconst. (11) 

2i 2t\ jjj J 

The only effect of scaling ^ is to change the unphysical additive constant in the action. 



2.1 Dynamics on a curved background M. 

To compute the large N effective potential and its minima, we need to evaluate the func¬ 
tional traces of equations (10) and (11) on Al. A natural way of doing this is by using zeta 
function regularisation. Recall that the zeta function of an elliptic operator A is defined by 


C(s) = TVA" 


Thus for the operator on M. we set 


In det 


-n +ml 


= — lim —Tr 
ds 




Differentiating this expression one obtains 


Tr [ ^ j = ^lirn^(sC(g + l))- (14) 

The main technical aspects of this work will involve the computation of the zeta function 
on a squashed three-sphere. The methods we use are similar to those in [24], although we 
are working in a nonperturbative framework. 
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3 0{N) model on a squashed three-sphere 


In this section we apply the above formalism with a squashed three-sphere as the back¬ 
ground. The metric of the squashed three-sphere of radius a is 

= T ^ > (15) 

where ai are the usual left-invariant SU{2) one forms 

(Ti -|- icJ 2 = {d9 + i sin 9d(j)) , = dip + cos ddcp , (16) 


and a is the squashing parameter which may take values in the range a G (—1, oo). In this 
parameterisation a = 0 is the round three-sphere. The other manifolds are also topologically 
spheres, but with a squashed fibration over 5^. The scalar curvature is found to be 


R = 


2(3 -h 4a) 


(17) 


0^(1 -|- a) 

Note that the scalar curvature changes sign and becomes negative when a < —3/4. We will 


also need the volume 


Vol(M) = 


2a^7r^ 


(18) 


(l + a)V2 • 

Although squashing the round sphere introduces an anisotropy, it preserves the homo¬ 
geneity of the space. Hence the spectrum of the scalar Laplacian on (15) is straightforward 
to calculate. As we saw above, this is what we need in order to compute the quantum 
correction to the self-energy which then enters the right hand side of the gap equation (10) 
of the large N theory on the squashed sphere. Following some minor manipulations [24], the 
zeta function (12) for the operator (—□ -|- m^)//i^ on the squashed sphere may be written 
as 

/(a/i)^^ 

\I2^ n(J _ 1 

1= 




tY “(^ — 1 — 2(^)2 -|- — 1 ]* 


(19) 


We may obtain a finite expression for ^(s) by analytic continuation of the sum (19). The 
methods are standard and are described in Appendix A. We find 


1 


d 


2s 


(is) = 


1 


-|- a 


m ^ 

2s 


2s 


+ 


4a' 


2s 


+ a^^e{a)H 


(3 -|- -|- a) 


~,2s 


r G{y)dy 

'o {1 + ay‘^y ' (1 +0)^70 e‘^^y + 1 


F{y)dy 




( 20 ) 


where the functions F,G and H are given in Appendix A and 0(a) is the Heaviside step 
function. We point out that although the appearance of the step function might suggest 



a discontinuity in the zeta function at a = 0, the function H and all its derivatives ac¬ 
tually vanish exponentially at a = 0, see equation (64) below. In fact the connection 
between the squashed sphere and thermal field theory is seen explicitly in the function H 
and dH(s)/ds\s=o (see Appendix A). The latter has a natural interpretation as the finite 
temperature free energy of a field theory on 5^ accompanied by certain parity reversals. 

For technical reasons in performing the analytic continuation, again see Appendix A, 
we have restricted ourselves to —8/9 < a < oo, omitting the range —1 < a < —8/9. The 
range we consider includes the value a = —3/4, where the Ricci scalar changes sign, the 
value corresponding to the round sphere, a = 0, and also arbitrarily large values of a that 
are relevant for the Hawking-Page type transition in the dual gravitational theory. From 
our results, the extrapolation to the range —1 < a < —8/9 seems straightforward, although 
the a ^ —1 limit itself is fairly singular and may contain interesting physics that deserves 
to be studied separately. 


3.1 Physics of the symmetric vacuum I: the mass gap 


In three dimensions, the coupling constant A has mass dimension one and thus the interac¬ 
tion is a relevant operator which generates a nontrivial renormalisation group flow from a 
free fixed point in the UV. High energies correspond to weak coupling, and at long wave¬ 
lengths the coupling becomes stronger. On the squashed sphere with scale size a a natural 
dimensionless coupling aX can be defined. We investigate the theory at different values of 
this dimensionless coupling, including arbitrarily strong and weak coupling regimes. Vary¬ 
ing this coupling naturally corresponds to varying the size of the compact space keeping A 
hxed. For any value of the coupling, the IR physics of the interacting large N theory is 
encapsulated in the solution to the mass gap equation (10). 

In three dimensions, zeta functions are regular at s = 1. Therefore the right hand side 
of expression (14) reduces to C(l)- The gap equation (10) then becomes 

A 


9 9 

m_ = m + 


rCm|(i). 


( 21 ) 


/iWolM 

To ensure that in the ultraviolet our theory flows to a free conformal fixed point, we set 
the bare mass term, m?, to be given by the conformal coupling to the background curvature 
R 


2 2 1 D ^ T 

= "^conf. = o^= ,0,. , ..A 


8 4a^(l-|-a) 


The gap equation may now be written as 


2 2 3 -|-4q; aA(l-|-a)^/^ Cm|(1) 

a m^ = ^ 


4(1 -I- a) 


27r2 


{afiy ■ 


( 22 ) 


(23) 
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The only dependence on a, the size of the squashed sphere, in this expression is contained 
in the dimensionless masses and couplings a\ and a?m^. 

It is clear that (21) always has solutions with > 0, essentially due to the 1/m^ 
behaviour in ((^ 2 ( 1 ) as ^ 0. In figure 1 we plot the solution for as a function of 

a for three fixed values of a\. The existence of a symmetric vacuum for all couplings a\ 
could only be seen because we have solved the theory nonperturbatively. At a perturbative 
level there is no symmetric phase for a < —3/4 and small oA. 



a 


Figure 1: a?m‘^ as a function of the squashing a with, from bottom to top on the left of 
the graph, a\ = 0.1, 10 and 1000. 

An interesting feature of figure 1 is that the value of on the round sphere, a = 0, 

is precisely the conformal value for a = 0, a?rn^ = 3/4, for all values of the coupling. This 
can be checked analytically using the expressions for the zeta functions in Appendix A. 
Thus the theory with tree level conformal mass term a?m? = 2)1 A posesses a self-consistent 
solution to the gap equation, for any non-zero coupling oA, precisely at = 3/4, which 

is again the conformal value! This can be understood as follows. On the critical 0{N) 
model at large N is obtained when vanishes. The round is conformally equivalent 
to M^. Thus it is natural that should take the conformal value on the round sphere. 
In fact this provides a check of our analytically continued expressions and will have a nice 
dual gravitational interpretation below. 

We now turn to the strong and weak coupling results for the theory which may also be 
interpreted as large and small volume regimes. 
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3.1.1 Weak coupling 

In the weak coupling regime we may solve for analytically. From the fact that C(l) ~ 
1/m^ as ^ 0 and from the gap equation (21) it follows that 

~ . for oA <C 1 if a < —, 

m^VolAl 4 

+ 0{X), for aA <C 1 if a > . (24) 

We can see that as aX 0+ the derivative of becomes discontinuous at a = —3/4. This 
is illustrated in figure 2. 



Figure 2: as a function of the squashing a < 0 and with aA ^ O'*'. 

The discontinuity suggests a phase transition at a = —3/4 in the extreme weak coupling 
limit aX 0'*'. Indeed ^ 0 in this limit. The vanishing of the “pion” mass, is 
reminiscent of symmetry breaking leading to massless Goldstone bosons. However, on 
a compact space there will be no symmetry breaking. Nevertheless, in the strict weak 
coupling limit there is a new phase that appears for a < —3/4. We know from (7) that at 
finite volume this must require (/>^ —> oo. Indeed we will see below that this phase transition 
is characterized by an order parameter, namely the condensate A < > which vanishes 

for a < —3/4 and is nonvanishing and of order N for a > —3/4. However, in most of this 
paper we will be interested in the case of finite coupling in which this phase transition is 
completely smoothed out. The nonperturbative smoothing of phase transitions may be a 
more general phenomenon of relevance in attempts to match perturbative phase transitions 
with strong coupling results [14, 15]. 
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3.1.2 Strong coupling 

For any finite value of the coupling the gap equation (10) seems difficult to solve analytically. 
In the limit of infinite coupling, the gap equation simply becomes 

Cm2 (1) = 0 as aX ^ oo . (25) 

Again, this equation does not appear easy to solve for general a but the roots can be 
determined numerically. However, it turns out to be analytically tractable at large a leading 
to a rather suggestive solution. Using the large a expression for the zeta function discussed 
in Appendix C we find that grows linearly with a 


The appearance of the golden mean in this expression is somewhat curious. Interpreting \foL 
as a temperature T, this equation gives us the high temperature behaviour of the solution 
to the gap equation, ~ T^. In fact, following the technique outlined in Appendix C one 
obtains exactly the same expression as above for the 0{N) model on x 5^ in the limit of 
high temperature, or shrinking S^. In this x limit, the expression (26) for the mass 
in terms of the golden mean was found previously in the 0{N) model in [25]. 

3.2 Physics of the symmetric vacuum II: the free energy 

Using the solution to the large N gap equation we compute the value of the action of the 
model as a function of a. We may do this keeping the volume of the squashed fixed 
while varying the squashing parameter. This has the effect of removing the logarithmic 
dependence of the action (11) on the volume of the squashed sphere. This action has a 
natural thermodynamical interpretation as the free energy times the inverse temperature. 
In figure 3 we plot the action of the symmetric vacuum (11) against a for strong coupling 
aX ^ oo. The qualitative form of the plot remains unchanged at finite values of the coupling, 
although it does change for weak coupling. 

The interesting features of this plot include the behaviour at large positive a, the be¬ 
haviour near a = — 1 and the presence of a maximum at a = 0. Based on numerical analysis, 
we are able to deduce that near a = —1, the action goes like I ~ —1/(1 + ot)^. For large 
positive OL the action scales linearly with a, which may be shown analytically. We will see 
below that all three features are also encountered in the dual gravitational theory. 



2 2 
a 
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TT^ 


In^ 


l + \/5 


as a 


oo. 
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Figure 3: The field theory action / as a function of the squashing a with aX ^ oo. 

We may find a closed expression for the free energy at large squashing a using (26) for 
strong coupling and (24) for weak coupling. The techniques discussed in Appendix C yield 

(27) 

and also remarkably 

(28) 

In these expressions Cr{s) denotes the Riemann zeta function. If one considers the squashed 
direction as a nontrivially fibred temperature, then the linear dependence on a is in¬ 
terpreted as a scaling with temperature. Writing F = IT, this gives the usual high 
temperature scaling of the free energy in a three dimensional spacetime. 

The factor of 4/5 difference between the strong and weak coupling regimes is reminiscent 
of the well known 3/4 factor that distinguishes the strong and weak coupling limits of the 
free energy of AA = 4 super Yang Mills theory in four dimensions [19]. In the present context 
we are able to calculate both limits nonperturbatively within field theory. The presence of 
Cr{^) in the free energies (27) and (28) is generic for conformal field theories on x 
the high temperature limit, see [25] and references therein. That paper also finds a factor 
of 4/5 in the free energy of the 0{N) model on x M^. 

4 The gravitational dual 

According to the original correspondence proposed by Klebanov and Polyakov [4] the grav¬ 
itational theory dual to the three dimensional critical 0{N) vector model is a higher spin 
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gauge theory in AdS 4 spacetime [26, 27, 28]. The details of this bulk theory, and hence the 
correspondence itself, have not yet been completely understood for AdS 4 spacetime whose 
boundary in global coordinates is the round three-sphere [29], For the 0{N) model on 
squashed spheres, the dual gravitational description should be a higher spin gauge theory 
on bulk geometries whose conformal boundaries are squashed three-spheres. We leave an 
in-depth study of the higher spin theory to future work. Instead, we will test the exact 
results from our field theory analysis against known semiclassical results regarding the as¬ 
sociated bulk geometries, without taking into account the effect of any higher spin degrees 
of freedom. We find a remarkably detailed qualitative agreement between the two pictures. 

In general there are two four dimensional Riemannian geometries with negative cosmo¬ 
logical constant that have the squashed three-sphere as conformal boundary. These are 
the AdS Taub-NUT and the AdS Taub-Bolt spacetimes. In Einstein semiclassical quantum 
gravity there is a first order phase transition between the two geometries at a critical value 
of the squashing parameter [5, 6]. This phase transition is simply the NUT charged version 
of the Hawking-Page transition for black holes in AdS [13]. More concretely, one finds 
that the AdS Taub-Bolt geometry only exists for a > 5 -|- 3\/3 10.2 and the quantum 

gravity phase transition itself occurs at acrit = 6 -|- 2\/l0 ss 12.3 [5, 6]. The authors of 
[5, 6] computed the difference of the actions associated to the AdS Taub-NUT and AdS 
Taub-Bolt spacetimes and inferred the existence of a transition at the point where the dif¬ 
ference vanishes. More interestingly from our point of view, in [7, 8] the individual actions 
for the AdS Taub-NUT and AdS Taub-Bolt geometries were obtained using a boundary 
countertem technique inspired by the prescription of [30]. We quote the result of [7, 8] for 
the bulk action as a function of the boundary squashing parameter 


Gtt (2a -|- 1) 
™ ^ ~GR (a+ 1)2 


(29) 


where G is Newton’s constant and R is the (negative) Ricci scalar for AdS Taub-NUT. Note 
that we use R to denote both the boundary and the bulk scalar curvatures. This solution 
exists for the entire range of a and dominates the partition function for a < 6 -|- 2\/T0. In 
particular near a = —1 the bulk action behaves as Jtn ~ ~1/(1 + ck)^ which is precisely 
what we found in the strongly coupled field theory. In addition Jtn has a maximum at 
a = 0, mirroring our field theory result of figure 3. 

For a > 6 -|- 2\/To the bulk action is that of AdS Taub-Bolt 


^TB = (1 + a) (^rub +‘^r{l + a) ^ , 


(30) 
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with 

and 

r = i(l + a)i/2 J'l ^ Vl-12(1+ «)-!+9(1+ a)-2) . (32) 

For large a, the Bolt action is negative and grows linearly 

Att 

Itb ~ a ^ oo . (33) 

This linear behaviour of the action translates into a scaling law with temperature, and 
hence a scaling of the free energy with temperature. This is essential in order to have 
an interpretation in terms of a three dimensional boundary field theory. Indeed, our field 
theory result (28) exhibits a negatively sloped linear dependence on a for large squashings. 
On the other hand, the AdS Taub-NUT solution vanishes for large a as /tn ~ 1/a. 



Figure 4: The bulk action / as a function of the squashing a. 

In figure 4 we have plotted the bulk action as a function of a. We have normalised 
the action so that its slope at large a agrees with the field theory result (28). There is 
an additive constant in the action which could be adjusted to make the peak values of the 
actions match. Comparing with figure 3, which was obtained from the strongly coupled field 
theory, we clearly see the qualitative similarities at small and large values of a. However, 
it is clear that there is no field theory analog of the Hawking-Page transition seen in the 
bulk. As we have pointed out earlier, we can presumably attribute this to the fact that the 
infinite massless higher spin degrees of freedom have not been taken into account in the 
bulk analysis. 
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We should emphasise that the qualitative behaviour of the field theory at very weak 
coupling is quite different to the strong coupling regime, particularly for the negative values 
of a where the boundary curvature becomes negative. At weak coupling we find that the 
action drops rapidly for a < —3/4 and scales as —1/(1 + which is faster than the 

gravity result which scales like —1/(1 + near a = —1. This is fairly straightforward to 
see following our treatment of the weak coupling regime. The strongly coupled held theory 
on the other hand, based on numerical evidence, appears to match the behaviour of the bulk 
theory near a = — 1. This suggests, as is usual in AdS/CFT dualities, that the strongly 
coupled held theory is dual to gravity on a weakly curved space. Below we will present 
another piece of evidence in favour of this identihcation which pertains to the behaviour of 
the held theory condensate (4> • <h) in the strongly coupled held theory and its implications 
for the gravity dual. 

Equating the bulk and boundary values for the effective action at large a, (28) and (33), 
suggests the tentative dictionary 


9Cij(3) N • 


(34) 


This formula passes the immediate test of giving weak curvatures when N is large. It might 
be compared with the well known result RG l/iV3/2 for the AdS^ X 5^ version of the 
correspondence. The appearance of Ci?(3) in (34) is curious and possibly tantalising given 
that Cr{3) also appears in a computation of a' corrections to the IIB effective supergravity 
action and hence the bulk free energy [31] for the AdS^ x version of the correspondence. In 
general we expect the coefficient in (34) to depend on the background. The same RG ~ 1/N 
relation was found for the Klebanov-Polyakov correspondence on AdS^ [29]. 


4.1 Negative curvature of the boundary 


Let us consider the physics associated with the sign change in the scalar curvature of the 
boundary at a = —3/4. From our comments above, the dual geometry to the squashed 
three-sphere in this regime is the AdS Taub-NUT spacetime which may be written as 


ds^ = 


1 4(12 , , 2^/ 2 , 2^ , + 2 

7271-272 ~r7- +r [l + ar ){ai+a 2 )+r — -^^3 


■4(1 


2 1 -b 


(35) 


where the cosmological constant is A = —3A;^ and the range of the radial coordinate is 
0 < r < 1. We have already argued above that the weakly curved gravitational description 
appears to be dual to the strongly coupled 0{N) model. Therefore we expect that the fact 
that a stable symmetric phase continues to exist for a < —3/4 at strong coupling should 
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translate into the bulk spacetime remaining stable in this regime. We have also seen that 
there is a phase transition at a = —3/4 for arbitrarily weak coupling. However, we expect 
this to be a strong curvature effect in the gravitational dual and thus beyond a classical 
gravitational computation. 

In string theory duals to field theories there are generically two types of instabilities 
that can arise when the conformal boundary has negative scalar curvature. For a summary 
of these, see [11]. The first is a nonperturbative instability due to the nucleation of BPS 
branes. The second is a perturbative instability due to negative modes of fluctuations about 
the solution. It is unclear whether the bulk higher spin gauge theory will contain objects 
analogous to BPS branes, so we will focus on the perturbative possibility. 

The perturbative stability of the AdS Taub-NUT spacetime against fluctuations of a 
scalar field was considered in [11]. That paper was primarily interested in scalar fields with a 
mass saturating the Breitenlohner-Freedman bound, as these arise from the compactification 
of eleven dimensional supergravity to four dimensions. The higher spin gauge theory that 
is conjectured to be dual to the 0{N) model has instead a conformally coupled scalar field, 
ip. It is simple to adapt the analysis of [11] to this case. We are interested in whether the 
bulk equation 

+(36) 
b 

has solutions with negative f3. This equation may be separated and converted into Schrodinger 
form 

=/3x , (37) 

using the transformations 

X = , dr = , (38) 

where g is the determinant of the metric (35). It is then easy to see that the potential H(r*), 
is everywhere positive for all values of the squashing a. Therefore /3 cannot be negative and 
the spacetime is perturbatively stable. 

One should also check stability against the other fields in the theory, such as metric 
fluctuations. However, stability of the bulk would be consistent with our field theory find¬ 
ing that the symmetric phase remains stable at finite coupling as we go to the regime of 
parameter space with a negatively curved boundary. 
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4.2 The condensate ($■$) 


Finally, we turn to an analysis of the condensate in field theory and its implications 

for bulk physics. We will also argue that the condensate can be used to construct an order 
parameter for the extreme weak coupling phase transition discussed in section 3.1.1. 

From the defining equations for the composite operator a and pion mass (3) and 
(9), we have that 

. (39) 

The Klebanov-Polyakov duality [4] then relates the vacuum expectation value ($•$) to the 
normalisable mode of the conformally coupled bulk scalar field (p. 

Let us firstly consider the strong coupling limit. We saw in equation (25) that rri^—m? = 
0(1) as a\ oo, that is, m\ remains finite as we take the strong coupling limit. Thus we 


have 


_($.$) ^ as aX^oo. 


(40) 


This implies that at strong coupling, aX ^ oo, the condensate vanishes. If we identify the 
strong coupling regime of the 0{N) model with the weakly curved gravitational dual, then 
this result translates into the (true) statement that the bulk field ip is not turned on in the 
AdS Taub-NUT and AdS Taub-Bolt geometries. 

The nonvanishing of the condensate at finite coupling then translates into the prediction 
that the bulk must allow for a nonvanishing value for (p in these cases. It is likely that these 
will just be new backgrounds with a nonvanishing p profile turned on self-consistently due to 
higher curvature corrections. Figure 5 shows a plot of ($•$) as a function of the squashing 
at two values of the coupling. 

One remarkable aspect of figure 5, which was anticipated already in figure 1, is that 
the condensate vanishes for the round sphere a = 0 for all values of the coupling. The 
dual implication of this fact is that when the bulk geometry is simply AdS, as opposed to 
AdS Taub-NUT or AdS Taub-Bolt, then it does not receive corrections at finite curvatures 
which turn on the scalar p. Such protection against corrections seems plausible, given the 
maximal symmetry of AdS, and mirrors known results for other realisations of the AdS/CFT 
correspondence, see for instance [32] for the AdS^ x version of the correspondence. 

At weak coupling the behaviour of the condensate (39) follows immediately from our 
previous results (24). In particular these imply that we may use the condensate as an order 
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Figure 5: The condensate as a function of the squashing a with a\ = 10 (top) 

and a\ = 1000 (bottom). 


parameter for the symmetry breaking phase transition occuring at a = — 3/4 as oA ^ O"*" 


lim 

aA^0+ 




0(1) for a < —3/4, 
0 for a > —3/4. 


(41) 


5 Some open questions 


Motivated by the duality conjectured by Klebanov and Polyakov [4], in this paper we have 
investigated and solved various aspects of the large N limit of the 0{N) vector model on a 
squashed three-sphere. Perhaps surprisingly we have found nontrivial agreement between 
our field theory results and the semiclassical gravitational physics of the proposed dual 
geometries: AdS Taub-NUT and AdS Taub-Bolt. 

There are various open ends to this work. The most obvious open question is the fate 
of the Hawking-Page transition in the bulk theory when massless higher spin degrees of 
freedom are taken into account. The field theory results suggest a smoothing out of the 
transition. 

Within field theory, one limit which we have not considered in detail is the a ^ — 1 limit. 
This is a rather singular limit, in which one dimension of the squashed sphere decompactifies. 
Nonetheless, perhaps there is interesting physics to be discovered. Our results near this limit 
have been mostly numerical. We have observed numerically that the strongly coupled field 
theory action in this regime scales in the same way as the action for the AdS Taub-NUT 
solution, with I ~ —1/(1 -|- a)^. An analytical solution of the field theory problem in this 
regime may shed further light upon the agreement between field theory and gravity. 
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Finally, it would be interesting to know if there are deeper reasons underlying certain 
amusing numerical coincidences we observed in the exact solution to the strongly coupled 
theory. Namely the appearance of the golden mean in the solution (26) to the gap equation 
at high temperature, or large squashing, and the appearance of 4/5 as the ratio of the weak 
and strong coupling limits of the free energy (28) at large squashing. These results are 
essentially due to the 0{N) model on x [25]. 
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Appendix A: Analytic continuation of the zeta function 


The zeta function we wish to analytically continue is 

oo l-l , 2s 


4iCW = EE 

^ 7—1 n 


la^ 


oo 1 — 1 

p+«('-1 - 29)2+Ajvif ^ E E^ 


(42) 


The method used has two steps, the first turns the q sum into an integral and the second 
deals with the I summation [24]. 

The Abel-Plana formula converts the q sum into integrals. To apply this formula we 
need to know the location of the branch points of f{q,l)- These are 

1 — 1 , i f P — 1 + a^rni ^ 


= 


a 


(43) 


We see that the sign of a is important for determining the position of the branch points. 
Let us firstly consider the case a < 0. In this case, the branch points are on the real axis 
and the Abel-Plana formula is 


q=0 




f{t, l)dt — 2i 


dt 


o27rt 


+ 1 


[f{l/2 + it,l) - f{1/2-it,l)] , 


(44) 


so long as the branch points q± are outside the range of integration (—1/2 ,1 — 1/2). For a 
given value of a this may be achieved by curtailing the range of 1. We will restrict attention 
to —8/9 < a which will allow us to use the previous expression (44) for all values I > 3. The 
remaining terms in the sum corresponding to I = 1,2 are then added to the final expression 
to obtain 




^ C(zz) = 4^ + 


4a 


2s 


2s 


m// (3 -|- -|- a)' 


+ a' 


2s 


t2s 


r G{y)dy 

0 (1 + ay^Y ' (1 + a)" Jq + 1 ’ 


F{y)dy 


+ 


(45) 
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with 


F{y) = E 


{p - A^y ’ 


r / i 

G{y) = ^ + iBy -e^y ~ {[I- wy -E^y 


where A, B,E are given by 


= 


B = 


1 — a?rrY 

1 + ay2 ’ 

1 + a’ 


2 ^ 1 - a mj 4ay^ 

1 + a (1 + aY 


Note that E'^ changes sign at an intermediate value of y. 

We now continue the sums in (46) using the Sommerfeld-Watson transformation. For 
the case of F{y) we have 

s if , z'^coiTTz 


= 5 X, *++# ■ 

where the contour Ci is shown in figure 6. By considering sufficiently large s we may rotate 
this contour into C 2 , also shown in figure 6, to obtain 

N if, z^cotvrz , 

F{y) = 2 Y ^ (™) 


-4 -3 -A -2 -1 


1 2 A 3 4 


Figure 6: Change of contours for the analytic continuation. 


This contour integral will require one more analytic continuation along the imaginary 


axis using 


coth nt = 1 + 


j27vt _ 
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We are interested in F{y) at s = 0,1 and in its derivative with respect to s at s = 0. One 
finds, after paying due attention to the branch cuts, the values 


F{y)\s=i = ^^cot7r^ + ^ 

p=i 

P{y)\s=o = - 5 - 


2 2 


j\2 _ p2 ^ 


(52) 

(53) 


The derivative at s = 0 may be shown to be 


A. 

ds 


Fiy) 


= 


s=0 


ittA^ 


+ A^ ln(l-e-2*^^) 


+ 


^ -2iTrA\ ^ 


TT 


T • / —2i7rA\ I T ' / —227rA\ 

■Li 2 (e ) + ^Li 3 (e ) 


+ ln|A^ 

p=i 


P 


(54) 


In this expression Li 2 and Lia denote the dilogarithm and trilogarithm functions, respec¬ 
tively. It is sometimes simpler to avoid the polylogarithm functions. An alternative expres¬ 
sion for when A is real is 

(■A ^ 2 


^F{y) 

ds 


s=0 


= —ttP f cot ntdt — In , 

■^0 27 r 2 ^ 


(55) 


whilst if A is imaginary we can write 


ds 


F{y) 


1^1 2ne , ittA^ Cij(3) 


5 = 0 J 0 


g27rt _ I 


dt + 


2tt‘^ 


+ ln|A^ -p^ 


(56) 


p=i 


In these expressions Cij('S) is the Riemann zeta function and P denotes the Cauchy principal 
value of the integral. 

Exactly the same manipulations are then performed for G{y). The contours Ci and 
C 2 are qualitatively as before, except that now the branch cuts emanate from the points 
2 ; = PE ± iB. The values are now 

2 


27r E sinh(27rR) — R sin(27rE) 


8Bp^ 


E cosh(27rE) — cos(27rE) |^2 _|_ ^^2 _ ^2 _j_ i2pB\‘^ 


G{y)l=o = 0 . 


(57) 

(58) 


The derivative is given by 
2 




s=0 


yy i2p In 

p=i 


{p -|- iB)"^ — E2 


— Fk 


B cosh 2ttB — t sinh 2nB — Be 
cosh 2nB — cosh 27rt 


— 27Tt 


{p - iBf - E2 

'^± -2iEln (^I - + -U 2 (l - 

- -Lia (1 - ) | + 2E In 


-dt 


cosh 2TrB — 1 


cosh 2TrB — cos 2TrE 


(59) 
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In this last expression, 1j\‘2{z) is the dilogarithm function. In evaluating these logarithms 
and dilogarithms, one should be careful to keep track of the phases of the complex numbers 
as E increases. In fact, the last two lines of (59) are given by 


—dvr 


B sin 2TTt — t sinh 2ttB 
cosh 2ttB — cos 2TTt 


dt , 


(60) 


when E is real and by 


—4Tr\E\B + dvr / 
Jo 


1^1 B cosh 2ttB — t sinh 2TrB — Be 


2'Kt 


-dt. 


(61) 


cosh 2TrB — cosh 27rt 

when E is imaginary. These integral expressions are easier to evaluate numerically than the 
dilogarithms. 

When a > 0 the branch points q±- are not on the real axis and = {I — l)/2 falls 
within the range of integration in (dd). Application of the Abel-Plana formula is as for 
a < 0 except that there is an extra contribution to the zeta function (d5) given by an 
integral along the branch cuts, which we take to run from q± to ±ioo, respectively [2d]. 
The extra term is found to be 


2^ da^"" . 

a H = — ——— sinTTS 


(da) 


E 

1=3 ■ 


I 


dy 


[1 + (y2 - p2)s ’ 


where 


p _ 1 + a^ml - 1 


a 


1/2 


(62) 


(63) 


Evaluating this expression for s = 1 one finds 


a^H 


Whilst the derivative at s = 0 is 



s=0 


^ 2aP [1 + (-l)'e2^^] ■ 


OO 

E 

1=3 


21 In 


l + (-l)'e 


I ^—2'kP 


(64) 


(65) 


The sums in the last two expressions are convergent and may be evaluated numerically. The 
last expression clearly displays the connection to a theory at finite temperature for a > 0. 
If we interpret y/a as a temperature, the above expression reduces to the finite temperature 
free energy associated to a field theory on a two sphere, along with certain parity reversal 
operations. 
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Appendix B: A check of the zeta function 

In [21] a series expansion in a is given for the zeta function on a squashed three-sphere 
for the value of the mass = R/6 (note that this is not the conformal value in three 
dimensions). The result is 

C{s) = (a/i)2^ [Ci?(2s - 1) + |sC/?(2s - 2) (66) 

+ s((/j(2s) -I- 2s(s -I- 1) — C_r(2s - 2) - — Ck(2s) —Cij(2s-|-2) ^ -|-. 

We have checked our expressions against this expansion and found very good agreement 
at each order in a. The agreement is rather nontrivial as it involves very precise cancellations 
between the terms involving F and G in the zeta function (45). It does not however provide 
a check of the H term (62) as this depends nonperturbatively on a. However we have derived 
this remaining term several times and furthermore it agrees precisely with the contribution 
calculated in [24]. 

Appendix C: The zeta function at large squashing 


The limit a ^ oo of the zeta function may be treated directly because one of the sums 
may be turned into an integral in this limit. The zeta function sum (19) may be rewritten 


cM—E E 


II -I-A-2-I- V 

k=-oo i£\k\+i+2n H— —j 

In the limit a oo, we may define x = and let 


E 

Zg|A:|+1+2N 


1/2 /-oo 


'(|fc| + l)/Ql/2 


This treatment is only allowed for the k = 0 terms if a^m^ja remains finite as a ^ oo. We 
will see that indeed this is the case at strong coupling, oA ^ oo. Doing the x integral, the 
zeta function becomes 


A/ N ^ _ 1 

4a^-i (s-1) 


as a ^ oo . 


To find the behaviour of at strong coupling and as a ^ oo we need to solve (25), 

that is C(l) = 0. We may evaluate C(l) by analytic continuation of (69) using the techniques 
of appendix A. The result is 


C( 1 ) = -^ 


a T^am-n 


2 a 


-k In I - 


2 _ g— 27ram7r/Q;^'^^^ 


as a ^ oo . 
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Setting this expression to zero one obtains equation (26) in text 


2 2 ® 1 2 
a ^ In 

TT^ 


l + \/5 


ctS CX 


oc, 


( 71 ) 


To compute the free energy in this limit, we need to evaluate C^(0)- By the usual methods 
we obtain 


C'(o) = 


a 


27r (am^ 


+ vr 


/ 


2{t^ — 


am-n I a 


c,2TTt 


- 1 


as a 


oo . 


(72) 


The integral may be expressed in terms of dilogarithms and trilogarithms. However, for the 
particular value of the mass given by (71) the expression simplifies dramatically to give the 
strong coupling result 

C^(0) = as a ^ oo , (73) 

Stt^ 

which immediately leads to equation (28) in the text. We have checked numerically that 
these strong coupling results actually match our analytically continued expressions involv¬ 
ing T, G and H in (45) in the large a regime. The matchings rely on delicate cancellations 
between the F, G and H terms and thus provide a nontrivial check of our analytic contin¬ 
uations. The weak coupling result (27) is obtained similarly, but in this case a^rn^ja 0 
as a ^ oo and one has to treat the /c = 0 sum separately. 
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